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ABSTRACT 


The  irreducible  binary  error  probability  of  an  FSK  system  with 
memory  operating  over  a  dispersive  channel  is  analyzed.  Binary  error 
probability  is  obtained  us  a  function  of  both  time-delay  and  Doppler 
spread  for  a  simple  HF  channel  model.  In  addition,  the  effectiveness 
of  correlated  diversity  receptions  in  combatting  channel  dispersion 
is  investigated. 
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I  INTRODUCTION 


This  memorandum  investigates  the  irreducible  error  probability  of 
an  FSK  communication  system  that  does  not  employ  energy  quenching  in  the 
detection  process.  In  contrast  to  the  incoherent  ma  t  c  h  ed  -  f' i  1 1  e  r  model 
(an  energy-quenching  device)  assumed  in  the  past,  1,2  our  model  more  closely 
approximates  existing  systems  that,  in  effect,  employ  narrowband  filters 
with  infinite  memory.  Thus  we  examine  system  performance  in  the  presence 
of  channel  dispersion  (at  infinite  s i gna l - t o- noi se  ratio)  when,  at  any 
given  time,  the  energy  stored  in  the  detection  filters  is  a  function  of 
the  entire  history  of  both  the  channel  variations  and  the  transmitted 
signal.  In  addition,  the  effectiveness  of  correlated  d ua 1  - d i v e r s i t y 
reception  in  combatting  channel  dispersion  is  investigated. 


FIG.  1  DUAL- DIVERSITY  FSK  COMMUNICATION  SYSTEM 


Figure  1  depicts  the  communiration  system  under  consideration.  Bold¬ 
face  notation  indicates  the  complex  envelope  of  the  corresponding  real 
narrowband  signal.  for  example,  the  real  signal 


1.  H.  F.  Daly,  "Analyaii  of  Multipath  Iffects  on  FSK  Error  probability  for  a  Simple  HF  Channel  Mode  1 ,  ** 
Heaearch  Memorandi  I,  Contract  SD-189,  SHI  Project  4554,  Stanford  Heaearch  Inatitute,  Menlo  Part, 
California  (Febriin,  1964). 

•!*  P.  A.  Bello  and  B,  iJ.  Nelin,  "The  Influence  of  Fading  Spectrum  on  the  Binary  Error  Probability  of 
Incoherent  and  Dt ' Cerent i a  1 1 y  Coherent  Matched  Filler  Here ivers, ”  I  HE  Trans,,  PGCS-10,  up.  lhO-168 
(June  1962). 


I 


*(  t) 


Rp  {«(()«■ 


is  the  on  •  l  lie  -  ii  i  i  t  lansmi  ss  i  on  of  an  KSK  system,  and 


J,(t)  Re  {Z,(<)e'2"/o'} 

is  tlie  real  signal  rereived  in  one  of  the  branches  of  the  du  a  1  -  d  i  v  e  r  s  i  t  y 
receiving  system.  The  diversity  channels  are  modeled  as  two  randomly 
time- varying  linear  filters  with,  in  general,  correlated  Gaussian  random- 
field  transfer  fund  ions,  //  j  (?,/)  and  H2  (  t  t  f )  •  3  If  X(/)  *s  the  Fourier 
transform  n  f  X(t).  then 


Z  t<f)  J  ,/>X</)e‘  277'fd/,  i  =  1,2 


(1  ) 


The  functions  h^(f)  and  (  t  )  are  the  complex  envelopes  of  the  impulse 
responses  of  the  narrowband  filters  centered  at  the  mark  and  space  fre¬ 
quencies,  /  -  i\  and  /()  +  \,  respectively, 

'Be'1 ' e'2”^{T- <  >  ,  t  >  0 


h  (  (  ) 


(2a) 


t  <  0 


h  ( t ) 


a 


t  >  0 
t  <  0 


(2b) 


After  narrowband  filtering  at  the  mark  and  space  frequencies,  each 
diversity  brunch  forms  the  difference  of  the  mark  and  space  square-  law 
detected  envelopes.  The  diversity  outputs  are  then  summed  and  sampled 
at  time  T .  yielding  the  decision  variable  Q . 


In  one  of  the  cases  to  be  investigated, 
tO 1  l  nw i n p  form: 


we  assume  X(t)  to  have  the 


X(  t  ) 


t  <  0 
0  <  t  <  T 
T  <  t 


(3) 


H.  I- .  I)al\,  4*  On  Modeling  t  hr  Time- Vary  inn  Kreque  nc  y-Se  1  ec  1 1  ve  Radio  Channel,"  Technical  Report  2— 
l*«iH  II,  (‘.nut  ran  DA  SC- WOHV* ,  SHI  project  4172,  Stanford  Research  Institute,  Menlo  Park, 

Cali  t<>]  uiu  (July  1 *M»4). 
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We  obtain  this  FSK  transmission  when  we  assume  mark  transmitted  in  the 
interval  (0,7’)  and  space  transmitted  for  all  f  outside  this  interval. 

At  time  t  -  T,  the  receiver  obtains  the  statistic 

Q  =■  |m,  (T)  | 2  -  |  st  <  r>  | 2  ♦  j  m2  ( r )  | 2  -  |s2m|2  (4) 

and  declares  mark  transmitted  in  the  interval  (0,7')  if  Q  0-  Thus, 
assuming  mark  and  space  to  be  equiprobable  messages,  the  probability  of 
error  P*  ,  for  X(f)  transmitted  is  the  probability  that  Q  is  less  than 
zero, 

p,  =  plQ  <  0]  .  (5) 


.  A. 
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II  CHARACTERISTIC  FUNCTION  OF  THE  DECISION  STATISTIC 


Decisions  are  based  on  the  statistic 


Q  -  Im, | 2  -  I s  ,  | 2  ♦  |m2 1 2  -  I s 2 1 2  .  (6) 

where  we  have  suppressed  the  argument  T  appearing  in  Eq .  (4).  The  sta¬ 
tistic  Q  is  a  quadratic  form  in  complex  Gaussian  random  variables.  Em¬ 
ploying  a  result  due  to  Turin,  4  we  obtain  an  expression  tor  the  characteristic 
function  of  Q . 

bet  V  be  the  vector  of  mark  and  space  filter  outputs. 


[m  jS  ,m2s  2] 


(7) 


and  (i  the  Hermitian  matrix 


then 


10  0  0 
0-100 
0  0  10 

0  0  0  -1 


Vt*GV 


(8) 


(9) 


is  a  Hermitian  quadratic  form  in  the  complex  Gaussian  variables,  flllf  Sj, 
IT10,  S  n .  Turin  showed  that  Q  has  a  characteristic  function  of  the  form 


n  (i  -  i  t\)~l 


(10) 


The  \ n  arr  the  eigenvalues  of  a  matrix  defined  as  follows:  Lee  Ky  be  the 

covariance  matrix  of  the  vector  Vt 


A  i 


SOL7**] 


(11) 


(J.  L.  Turin,  "Thr  Characteristic  Function  of  Hermitian  Quadratic  Forma  in  Complex  Normal  Variables, 
fitoaefrtfcu,  Vol,  47,  pp.  199-201  (June  i960). 
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i 


then  the  Kn  in  Kq .  (10)  are  the  eigenvalues  of  the  matrix 


KVG  -  E[VVt*G] 


m  jin* 

-m,s* 

m  j  m  * 

-m  ,s 

s,«*. 

-s.s* 

S.m2 

-$  jS 

m2ni* 

'm2S* 

m2m* 

-m2s 

s2m; 

-S2S’ 

s2m*2 

-s2s 

(12) 


In  this  memorandum  we  evaluate  the  eigenvalues  of  the  matrix  K  yG, 
invert  the  characteristic  function  <Pq(  t ) *  and  obtain  an  expression  for 
the  probability  of  error  P 
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Ill  DIVERSITY  CHANNEL  MOOELS 


Two  correlated  Gaussian  random-field  transfer  functions, 
and  Hn{t,f),  serve  to  model  the  diversity  channels.  We  assume  that  the 
random  fields  are  homogeneous  and  identically  distributed, 

£[«*(  t  *  a,f  ♦/?)]  -  ♦  <*,/  ♦£)]  =  fl„(a,/3)  (13a) 

and  that  the  cross  correlation  between  the  channels  is  simply  given  by 

E[H\(t  ,f)H2(t  *a,f  *fi)]  =  pRh(ol,/3)  ,  (13b) 

whe  re 

o  <  ip  I  <  1  . 

Observe  that  Eq.  (13b)  can  be  interpreted  in  terms  of  a  three-dimensional 
random  field  (for  example,  we  may  consider  diversity-antenna  separation, 
d,  as  the  third  dimension)  in  which  the  t h i rd- d ime n s i ona 1  dependence  can 
be  factored  from  the  three-dimensional  autocorrelation  function  ( d ,  0£ ,  /3 ) 

ftw<<i,a,/3)  =  C(d)R„(a,t3)  =  pHH(  a,/3)  .  (14) 

The  diversity  correlation  coefficient  p  as  formulated  in  Eq.  (13b) 
is  defined  in  terms  of  complex  envelopes  and  can  itself  be  complex.  In 
practice,  one  usually  deals  with  the  correlation  between  real  envelopes. 
In  what  follows  we  relate  the  magnitude  of  p  to  an  easily  measurable 
quantity,  the  correlation  coefficient  between  the  squared  envelopes  at 
zero  time  lag  and  zero  frequency  shift: 

Co  v  [  I  // ,  |  2,  |tf2|  2] 

y  -  -  (15) 

Var^tlff,  I  2] Var**[|/f2|  2] 

where  the  I  and  f  dependence  has  been  suppressed  in  the  H 
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Cov  [|tf,  |  2,  \h2\ 2)  -  EI\h,\2\H2\2]  -  £[|w,l  2]£[|//2I  2] 

£ [  I W,  I  2|//2I  2]  =  ,tf2//*2] 

=  £[//*//,  ]£[ff2«2]  ♦  £ttf*tf2]  £[//,//.,]  ♦  £ t//,// 2 1 

=  rt2(0,0)  ♦  0  +  |p|  2fl2(0.0) 

=  (1  *  l/o|  2)fl2(0,0)  . 


Cov  [  |  W,  I  2,  I«2|  2] 
Var  [  |//2  |  2] 

£[|W,  I  4] 


Var  [Ilf,  I2] 


(1  ♦  Ipl  2)flj|(0,0)  -  fl2(0,0)  ■  lpl2fl£(0.0)  . 

Va r  [  I  //,  I  2 ]  «  £f  I W,  |  4]  -  £l  I  //,  I  2]  . 

£[//,//*]£[//,//;]  *  £[«,//,]£[«*//*]  ♦  £[//,//;]£[//,//;] 

2fl2(0,0)  . 

2fl2(0,0)  -  fl2(0,0)  =  fl2(0.0)  . 


From  Eq .  (15),  we  obtain 


P  2«£<0.0) 

y  ■  - - -  -  \p\  2 

R2„(0,0) 


(16) 


To  simplify  the  evaluation  and  interpretation  of  the  FSK  error- 
probability  performance,  we  restrict  our  attention  to  the  case  in  which 
P  is  a  real  quantity.  In  this  case,  the  error  probability  will  depend 
on  p  ^ ,  the  squa red- e nve lope ,  diversity  correlation  coefficient. 


✓ 
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IV  FILTER  OUTPUT  COVARIANCES 


The  elements  ol  the  matrix  K  VG  (omitting  sign)  are  the  elements  of 
the  matrix  Kv ,  the  covariance  matrix  associated  with  the  filter  outputs 
In  lieu  of  the  prior  noLation,  III  ^  IH2,  S2>  it  is  convenient  at  this 
point  to  express  the  four  filter  outputs  as  follows, 


V  (r) 

t 


f  7 

Z,(  ( )h  r(T 

J—  OD 


t  )dt 


1,2 

m  t  s 


(17) 


where  the  index  i  indicates  the  diversity  branch  and  r  the  mark  or  space 
channel.  All  terms  of  the  covariance  matrix  can  be  written  in  the  form 

r*  fT 


£[Z*( ( l)Zj(t2)]h*  (T  -  t  ,)h 


,(T  - 


J—  an  J—  cd 


(18) 


Equations  (1)  and  (13)  imply  that 


£lz*  ( r,  nj  l  tj)I  -  t.  t]j  Jt't  f 


T)K(  (,  - 


i  2wA  <  t . 


t i dhdr 


i\9) 


w  In-  r  t* 

U  M  .  - 

J 

4  ■  < 

\p  *  r  ‘  * 

j 

a  n  d 

f  f  “  i  2^r  (Aa-r/3  ) 

( at/3)d<xd/3 

S,iK,r>  -  jj  c  rh 

is  the  channel  scattering  function. 
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The  symmetric  two-path  case 


S„(\,T) 


fl„(0,0)  f 

-  <5(\  -  K'  )S(r 

2  l 


T  ') 


•f 


b (K  ♦  K' )b{r  ♦ 


(20) 


is  particularly  useful  in  analyzing  systems  operating  over  HF  channels. 
The  channel  is  modeled  as  two  discrete  paths  spaced  2 T*  in  time  delay 
and  2 Kl  in  Doppler  shift.  Let 


A 


x  r 


(K,t) 


r 

X(t  -  T)h r(T  -  t )e ,2wA  'dt  , 

—  00 


m ,  s  , 


(21) 


then,  lor  the  symmetric  two-path  case, 

Vo.o)  f  1 

£[V*(r,)r  ( P,)J  -  ktj  — —  6:  <V,t'  >4  (a\t* )  (-V,-r*  M  <-V,-r' ) 


(22) 
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V  THE  EIGENVALUES  OF  K„G 


Le  L 


1  r  2 


E[Kr(ri)v.(r2)1 


lor  r (  =  m ,  s  ,  r  2  =  « , s  and  i  s  1,2  ,  and 

-c 


i  a  s  s 


then  it  follows  from  Eq .  (22)  that  the  matrix  KyG  can  be  written 

r  n 

C  pc 

i^c  = 

pc  c 

We  seek  solutions  (\  and  K)  to  the  equation 


KVGY  * 


By  setting 


and  using  Eq .  (24),  Eq .  (2  5)  becomes 


PC 

F, 

=  K 

i'i 

1 — 

rj 

c 

^2 

Let  the  matrix  C  have  eigenvectors  X%  and  eigenvalues  co .  ,  i  = 


t  he  n 


cx 


<^>tXl  .  i  *  1,2 


(23a) 

(23b) 

(24) 

(25) 


(26) 

1,2. 

(27) 


to 


In  Eq.  (26)  let  K,  *  Y2  =  Xt  , 


"(1  ♦  p)CX , 

=  K 

(1  +  p)CXl 

(1  +  p)u>iXl 

~X,~ 

=  K 

(1  +  p)u  x, 

_X,_ 

thus,  if  C  has  eigenvalues  ojj  and  u>2  ,  then  KyG  has  eigenvalues  (1  +  p)a>j 
and  (1  +  p)oo2.  In  Eq .  (26)  let  Y{  -  -Y2  =  X  t 


(1  -  p)CX, 

X  ’ 

1 

=  X 

-(1  -  p)CX , 

-x>_ 

(1  -  p)cotXt 

X 

l 

=  K 

-U  -  p)o>,Xt 

-X 

1 

thus,  if  C  has  eigenvalues  co j  and  o>2  ,  then  K  yG  has  eigenvalues  (1  -  p)co  j 
and  (1  -  p)  cl>2  .  Let  the  four  eigenvalues  of  KyG  be  denoted  by  i  *  1,2, 

3,4  and  the  two  eigenvalues  of  C  by  oj  ,  i  •  1,2,  then: 


\  * 

(1  ♦  p)co  J 

(1  -  /D)w, 

s  - 

(1  ♦  p)c^2 

> 

11 

{  1  -  p)u>2 

(28) 

To  evaluate  6Uj  and  u>2 ,  we  fi 

nd  the  roots  of  the  quadratic 

equation, 

lc 

-  col  |  *  0 

^  -  (c..  ■  C.s)aJ 

♦|c  |2  -  c  c  .  =  0  , 

(29) 

“l  <c..  •  c,.  )  1 

2  "2 

|^c..  +  c,.  )2  '  ^ 

(30) 

11 


Out*  nni  verify  that  thi*  above  roots  are  real  and  that  a>j 
,,  negative.  The  eigenvalues  and  are  real  and  posi 

eigenvalues  V  ^  and  A  ure  real  and  negative. 


i s  pos i t i v  e  and 
t ive ,  and  the 
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VI  ERROR  PROBABILITY  EVALUATION 


The  probability  density  of  the  decision  statistic  Q  is  found  by 
inverting  the  characteristic  function, 


V<7> 


'q<PQU  )dt 


fa <<?> 


e  *  q<t>Q  ("  i  s  )ds 


l  r+  |C 

277!  J 


'«  n  (1  -  \  s)~lds 

n  ~  l 


For  (7  <  0 ,  the  integrand  in  Eq  .  (31)  goes  to  zero  along  an  infinite  semi¬ 

circle  in  the  left-hand  plane;  for  q  >  0,  the  integrand  goes  to  zero  along 
an  infinite  semicircle  in  the  right-half  plane.  Thus,  letting  rn  denote 
the  residue  of  the  integrand  at  the  pole  s  =  1  /K  and  evaluating  the 

integral  in  Eq  .  (31)  by  residues,  we  obtain: 


,  for  q  >  0 


r3  ♦  r4  ,  for  q  <  0 


n  (\n  -  Kk)~l 

n  t  k  n 


(The  eigenvalues  and  \2  are  positive,  and  the  eigenvalues  A3  and 
are  negative. )  Thus,  given  that  mark  is  transmitted  in  the  interval 
(0,T),  the  probability  of  error  can  be  written 


13 


P.  PiQ  <  0]  -  /*  fQlq)dq  , 

P  II  (\  -  -  \J  n  (A.  .  (34) 

nN  nr  4 


The  expression  for  the  probability  of  error  can  be  put  in  terms  of 
the  elements  of  the  matrix  C  and  the  parameter  p.  Let 


ai,  - 


< < cma  *  c,,)2  *  4Ic«,I2> 


K 


(35a) 

(35b) 


then 


(Y  -  X)2  [(1  *  p2)Y  *  2X] 
4x7P  -  p*Y2] 


(36) 
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VII  THE  FILTER-SIGNAL  CROSS-AMBIGUITY  FUNCTIONS 


The  functions, 

4,„(\,t)  =  fTX(t  -  T)hm(T  -  t)e‘2nA'dt  (37a) 

-  CD 

=  /  X(t  -  r)h,(T  -  t)e,2nA,dt  ,  (37b) 

-  CD 

appearing  in  Eq .  (22)  represent  the  outputs  of  the  mark  and  space  filters 

at  time  T  when  the  channel  introduces  a  time  delay  r  and  a  Doppler  shift 
In  this  section  we  document  these  functions  for  two  different  trans¬ 
missions  and  the  narrowband  filters  assumed  in  Eq .  (2).  Since,  for  the 

filters  assumed  in  Eq .  (2), 


^,(X,r)  =  Asm{K  -  2A,t) 


it  is  sufficient  to  display  4  (\,r) 

Let 


x(0 


t  2T7A( 

"  «  2  77  A  t 


i  2ttA<  t-2D 


,  t  <  0 

.  0  <  I  1  T  , 

,  T  £  t 


then 


/t,.U,T) 


ABe~aT 

a  7  i277(  A  ♦  2A) 


exp {  [a 


i2t7(A  ♦  A)]t> 


(38) 


AB  exp{-aT  ♦  i2t7At> 
a  +  i2^ 


[exp{(a  ♦  i2rrK)  T } 


exp ( ( a  ♦  i  2vK)r] ]  , 


for 


0  <  r  <  T 


(39a) 
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from  Eqs.  (39a)  and  (39b). 

The  error  probability  behavior  is  independent  of  any  normalization 
of  the  function  4^h(X,t).  It  is  also  interesting  to  observe  that 
4  (A.,t)  depends  only  on  the  following  normalized  parameters: 

OlT,  a  normalized  filter  time  constant, 

A  T,  a  normalized  frequency  shift, 

XT,  a  normalized  Doppler  shift, 
r/T,  a  normalized  time  delay. 
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VIII  CONCLUSIONS 


Figures  2  through  9  summarize  the  asymptotic  or  irreducible  error- 
probability  behavior  of  a  du a  1 - d i ver s i ty  FSK  system  with  memory.  All 
figures  pertain  to  the  case  in  which  the  mark  and  space  frequencies  are 
separated  by  the  inverse  of  s i gn a 1 i ng- e 1 emen t  duration  ( 2A T  =  1).  The 

symmetric  two -path  case,  described  in  Eq.  (20),  is  assumed.  The  differ- 
erence  in  time  delay  between  the  two  paths,  divided  by  T,  is  referred  to 
as  the  time-delay  spread  2crr.  Similarly,  the  difference  in  Doppler  shift 
between  the  two  paths,  multiplied  by  T,  is  referred  to  as  the  Doppler 
spread  2  cr^.  The  transmission  described  by  Eq .  (38)  is  referred  to  as  the 
infinite  space  past  case  and  the  transmission  described  by  Eq.  (40)  is 
referred  to  as  the  one  element  space  past  case. 

It  is  interesting  to  observe  that  the  time  constant,  OLT  =  2,  yields 
a  significant  improvement  in  performance  over  that  obtained  for  OLT  -  1. 

Also  of  interest  is  the  fact  that  diversity  is  effective  in  combatting 
channel  dispersion  for  correlations  as  high  as  0.8. 
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FIG.  2  ASYMPTOTIC  BINARY  ERROR  RATE  vs.  2rrr  INFINITE  SPACE  PAST,  a  T 
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FIG.  3  ASYMPTOTIC  BINARY  ERROR  RATE  vs.  2<ta  INFINITE  SPACE  PAST,  «T  =  1  AND  „T 
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FIG.  4  ASYMPTOTIC  BINARY  ERROR  RATE  vs.  2<rr  ONE  ELEMENT  SPACE  PAST,  «T  =  1  AND  «T 
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